Inclusive scattering data on light nuclei as a tool for the extraction of G 
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We demonstrate that refinements in the analysis of inclusive scattering data on light nuclei enable 
the extraction of, generally accurate, values of the neutron magnetic form factor GJf(Q 2 )- In par- 
ticular, a recent parametrization of ep inclusive resonance excitation enables a reliable calculation 
of the inelastic background, and as a consequence a separation of quasi-elastic and inelastic contri- 
butions. A far larger number of data points than previously considered is now available for analysis 
and enables a more reliable extraction of G^.j from cross section and Rt data on D and He. The 
achieved accuracy appears mainly limited by the present uncertainties in the knowledge of proton 
form factors and by the accuracy of the data. 



INTRODUCTION. 



In a previous report we discussed the feasibility to extract the neutron magnetic form factor (FF) G , ^ I (Q 2 ) from 
inclusive electron scattering data on D and 4 He [j| at small and moderate Q 2 . Using a simple model for inclusive res- 
onance excitation, we established, that in the immediate vicinity of the quasi-elastic peak (QEP), the QE components 
dominate the inelastic background, which was subsequently neglected. That procedure severely restricts the number 
of data points around the QEP available for the extraction of a n — G^/ ' \x n Gd (Mn an d Gd{Q 2 ) = [1 + Q 2 /0.71]~ 2 
are the static magnetic moment of the neutron, and the standard dipole form factor). 
For several reasons we propose to extend the above analysis: 

1) Availability of an unpublished set of low-Q 2 quasi-elastic (QE) data on 4 He, 2], which are contained in a PhD 
Thesis of J-P. Chen Q and which cover a nearly continuous range in Q 2 . Moreover, contrary to the older NE3 data 

most of the NE9 data sets reach, or extend into the resonance region. 

2) After completion of Ref. {I] we have been informed of recent high-quality inclusive resonance excitation data on 
a proton, which supersede parameterizations of older SLAC p| and of more recent JLab data @ (see also Ref. 0)- 
In addition a model has to be devised in order to obtain the required neutron structure functions (SF) FJ} (k = 1, 2) 

which appears in the desired F? , the averaged p, n SF. The latter is related to the nuclear SF F k by means of 
f^ N ' A , the SF of a fictitious nucleus composed of point-nucleons, which has to be computed. 

The above F£ may be separated into nucleon-elastic (NE) and nucleon-inelastic (NI) components F k ' , F k ' , 
which correspond to processes where, after absorption of a virtual photon, the nucleon is not, or gets excited. The 
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above mentioned relation generates a similar division of the nuclear SF F k = F^ NI + F£' . 
In each data set we consider two regions: 



, — i 

X. 

A) The high energy-loss region far from the QEP, which is dominated by the nuclear NI background: Tails of QE 
contributions barely depend on the precise value of G^j and, the overwhelmingly inelastic computed cross sections 
can thus be compared against data. 

B) With decreasing energy loss one reaches the inelastic and elastic sides of the QEP. There we focus on the 
difference between data and the NI components. To those we apply a criterion which, when fulfilled, identifies the 
above differences as F k A ' NE , which is related to the nucleon FFs. The QE region will be shown to be the most sensitive 
one for variations in cv„, and is hence the area of prime interest. 

The above program is hampered by complicating circumstances. Foremost is the discrepancy 0, ITo| between the 
ratio, extracted from Rosenbluth-separated elastic ep cross sections and from e(p,p)e polarization transfer 
data [jj], At this time it is not obvious which p FFs should be used. It seems that two- photon exchange 

contributions [l3L[l3| confirm the E/M ratio for the proton, provided by the polarization transfer data. Unfortunately 
this does not directly reflect on the FFs themselves, because the ep data have as yet not been corrected for those 
contributions, prior to an extraction. Before resolution of the above issue, a choice for the required input will have to 
be made. Arrington recommends the use of the experimentally simplest, non-separated ep cross section data. 

The present note is organized as follows. We start with total inclusive cross sections, but shall also re-analyze the 
components for the absorption of virtual transverse photons ex Rt- After discussing the required input, we analyze 
all good-quality data on D, 4 He. 

We first concentrate on the region below and around the QEP, where NE components generally dominate. This 
enables one to reliably extract a n from any given data set. We no more neglect NI, as had previously been done 0, 
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and instead subtract those from the data in order to isolate the NE components. We also show that the computed 
cross sections agree very well with the data in and beyond the resonance region, which are completely dominated by 
NI components. Between the inelastic slope of the QE region and the elastic wing of the resonance region, the NI 
components falls short of their predictions and we discuss a possible origin. 

In the discussion we set limits to the accuracy of the results, caused by occasionally data of insufficient quality, the 
ambiguity of the proton FFs, and uncertainty in the NI components. Without substantial improvements in each of 
the above items, we do not envisage the possibility to substantially improve the results of the model. 



II. QUASI-ELASTIC INCLUSIVE SCATTERING. 

Consider the cross section per nucleon for inclusive scattering over an angle 9 of unpolarized electrons with energy 



E 



d 2 <7 A (E;6,v) 
dfl dv 



a M (E;9,is) 



^F A ( X) Q 2 ) + ^F A (x,Q 2 )tg 2 (e/2) 



(2.1) 



with CM, the Mott cross section. Above, F A 2 (x,Q 2 ) are the nuclear SF, depending on the squared 4-momentum 



transfer q 2 = -Q 2 = —(\q\ 
< x < A. 

In an alternative representation 



v 2 ) and the Bjorken variable x — Q 2 /2Mv, with M the nucleon mass. Its range is 
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dfl dv 



= a M [Wf + 2W{-tg?(e/2)\ 
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(2.2) 
(2.3) 



where e -1 = 1 + 2(|q| 2 /Q 2 )tan 2 (6 l /2). The SF W A , W A <|2.3j) relate to the absorption of longitudinal and transverse 
photons. 

In the sequel we shall use a relation between SF F k ' A for nucleons (N — p, n) and a nucleus which, for isospin 
1 = targets of our interest, reads 0] 



F k A (x, Q 2 ) = ( A -^f PN ' A (z, Q 2 ) ]T C kl (z, Q 2 

J X z , 



(2.4) 



The link between F N = (F p + F n )/2 and F A is provided by f PN ' A , the SF of a fictitious target A, which is composed 
of point- nucleons. Eq. 1)2. 4|1 holds in the Bjorken limit Q 2 — > oo, as well as in the Plane Wave Impulse Approximation 
(PWIA). 

For quite some time we have considered Eq. (|2.4|) for finite Q 2 in an alternative, non-perturbative approach |15| . 
based on a covariant generalization |l6( of the non-relativistic theory of Gersch-Rodriguez- Smith (GRS) ml We 
considered the above as a conjecture, and its apparent validity for Q 2 > Qq ps 2.5 GeV 2 as an empirical fact |18|. 

Only recently did we come across work by West and Jaffe who more than 20 years ago proved Eq. (|2.4|l in the 
PWIA, using either a parton model or pQCD [T^.l20l |. It is actually possible to generalize their proof by adding Final 
State Interactions (FSI) to the PWIA, reaching the Distorted Wave Impule Approximation. The same holds for the 
inclusion of FSI in the GRS version [21| and the proof formally recovers Eq. I|2.4|l . The intriguing difference lies 
in the interpretation: in the effective hadronic description one uses typical nuclear concepts, as are nuclear density 
matrices, effective iViV scattering amplitudes, etc. Those are of course foreign concepts in QCD. Similar remarks hold 
for recently discussed effective nuclear parton distribution functions [22|. 

We return to Eq. 1)2. 4fl . which includes the effect of mixing of the nucleon SF embodied in the coefficients Cm 
[23Ll24]| . In both the PWIA and the GRS approach, C\\ = 1, C12 = 0, while C21 is negligibly small. For a discussion 
of an approximate fashion to compute C22 in the GRS, we refer to Appendix A of Ref. Since the approximately 
calculated deviation of C22 from 1 does not decrease fast enough with Q 2 , we use C*22(Q 2 ) — > 1 for Q 2 > 3.5 GeV 2 . 

Finally we remark that Eq. I|2.4|l relates to nucleons as the dominant source of partons. This is the case for x > 0.20 
[25| and thus certainly for the range on which we focus 0.4 < x < 1.2, which comprises the QEP. 

Next we recall the NE and NI parts F k N ' NE , F^ 
virtual photon onaiV 7* + N — > N and to inelastic absorption 7* + N 



pN.Ni Q £ nuc j eon gp^ w hich correspond to the elastic absorption of a 



(hadrons, partons). Elastic components 
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for a nucleon vanish except for x = 1 and contain the standard combinations of the electro- magnetic FF G E M (Q 2 ) 
Denoting the average of their squares by [G N ] 2 = [(G p ) 2 + (G") 2 ]/2, one has (r? = Q 2 /(4M 2 )) 



F?>^(x,Q 2 ) = -6(l-x){GZ(Q 2 )} 



F, 



N,NE 



(x,Q 2 ) = 6(1 -x) 



m(Q 2 )} 2 

iG^mr+rtGUQ 2 )} 2 } 



(2.5) 
(2.6) 



It is a trivial matter to express the corresponding NE, NI parts of nuclear SF using the link 12.4(1 in its region 
of validity x > 0.2. Thus for the nuclear NE (QE) parts 

\x,Q 2 ) = fPN,A< f' Q2) [G N M (Q 2 )] 2 ] 



pA,NE, , ,2 , 



F 2 A - NE (x,Q 2 ) = xf PN > A (x,Q 2 )C*(x,Q 2 ) 



g n e {Q 2 )\ 2 + ii[g n m {Q 2 )] 2 

A,NE 



1+7? 



(2.7) 
(2.8) 



Since y N - A for light nuclei is sharply peaked around x ~ 1, the same holds for F k 
corresponding to Eqs. (2.5), (2.6) one has 



N.NE 
T 

N,NE 



f/nrJV 12 

( Xl Q 2 ) = 6(1-x)^ m1 



(x,Q 2 ) = 6(1 -x) 



M 

(l+r/)[Gff] 2 
2Mn 



For the L, T components 

(2.9) 
(2.10) 



Clearly the NE L,T components for the N separate magnetic and electric FF. However, the generalization of the 
above to composite targets depends on the model which relates NE parts of the nucleon and target SF. In the GRS 
approach one finds from Eqs. (2.4), (2.9) and (2.10) 
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(2.11) 



(2.12) 



The NE part of the transverse nuclear SF still contains only magnetic FFs. However, for the model defined by 



+ {(C 22 -l)(x 2 +n)+x 2 -l\[G N M ] 2 



Eq. (|2.4|) the longitudinal nuclear partner R 



A,NE 



Eq. (|2.12|) is generally a combination of E, M FFs, except at 



the position x — 1 of the unshiftcd QEP, and then only for unmixed Fj? in Eq. (2.4), i.e. with C22 = 1- The 
above contrasts with the PWIA, where no such mixture occurs in R^' NE . Only in a very limited number of inclusive 
scattering experiments have Rosenbluth L,T separations been performed. As a consequence most of our efforts are 
concentrated on the more involved total inclusive data. 

When focusing on FFs, one has to isolate in the data the NE parts (2.7), (2.8), which contain those FFs. Such a 
procedure obviously requires accurate knowledge of the nuclear NI background in the QE region and which on 

the adjacent inelastic side of the QEP x < xo(Q 2 ) < 1, is dominated by inclusive resonance excitation. 

Just as the NE components in F, ' E produce in nuclei the corresponding F k ' NE , which is centered around the 
QE peak, also A^-resonances in F k ' NI may be reflected in nuclei SF F k ' NI . In particular for A < 4, for which f PN < A 
is sharply peaked in x, both F^' 1 ^ 1 and F k ' NI have maxima at about the same x. We recall, that the observed 
structures in F k ' NI for Q 2 < 3 GeV 2 and v > i^qep, are not genuine target resonances but nucleon resonances, 
modified by the nuclear medium [2(|. We refer to those structures as 'pseudo-resonances' |27j . 

We conclude this section returning to the validity of Eq. (|2.4(l . The simple forms (12.7(1 and 1(2.8(1 for the nuclear 
NE parts hold quite generally, even below the estimated Q\ ~ 2.5 GeV 2 . The same may be assumed for the NI part, 
due to the inclusive excitation of very narrow resonances. With some hesitation we shall therefore make applications 
even for Q 2 = 0.5 4- 1 GeV 2 , but one should be prepared to encounter less good fits for those, than for larger Q 2 . 
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III. INPUT. 



We review major theoretical input elements: 

1) Density matrices for the target nuclei, diagonal in all coordinates except one. Those relate to ground state 
wave functions and have for the lightest nuclei been calculated with great precision [2^, |2j| . For heavier targets one 
has to invoke approximations, for instance by interpolation between diagonal density matrices and special limiting 
situations [l7l.l30| . Wishing to avoid theoretical uncertainties, we do not incorporate in our analysis data for targets 
with A > 12. 

2) NN dynamics for Final State Interactions (FSI) [30|. which enters the calculation of f PN ' A . 

3) Ff' : We considered various representations for Ff , all having explicit resonance and background components, 
namely Stuart etal p| and more recent ones 'christyl' and 'christy2' Q, based on Rosenbluth-separated cross sections 
(cf. also Q). The second version is claimed to be of somewhat better quality. However, the argument x/z oi the 
latter in Eq. I|2.4() varies for fixed Q 2 and occasionally crosses stated regions of validity. The ensuing inconsistencies 
appear more severe for christy2 than for christyl and we thus prefer the latter. Since the parameterizations are of 
relatively poorest quality for small Q 2 , one should expected correspondingly inferior results. 

For large Q 2 one avoids uncertainties and even inconsistencies by choosing a fixed Q 2 = 3.5 GeV 2 , beyond which we 
switch to a parametrization of F p , averaged over resonances |3l|. For Ff at high Q 2 we employ those of Bodek- Ritchie 
[32|. The chosen procedure is in line with quark- hadron duality, which predicts similar outcome for Fj^ , both globally, 
when averaged over the entire resonance region |3^| , and locally for isolated resonances 0, |35| . 



4)F, 



n:NI 



We use a procedure presented in Ref. |8|, based on the ratio 

C(x,Q 2 ) = F 7l (x,Q 2 )/F p (x,Q 2 ) 



J2d k (Q 2 )(l-xY 



k=0 



with the coefficients dk{Q 2 ) to be determined by information on C for 3 selected points: 

i) C(0, Q 2 ) — 1, required to obtain a finite Gottfried sum 

ii) Use of the primitive choice F n = 2F D — F p , which is accurate for x < 0.3 and is exploited for the value x 

iii) Information from the elastic end-point x — 1 , where C is determined by static FF 



C(1,Q 2 ) 



~] 2 



G n E (Q 2 ) 



GIAQ 2 ) 



G P E (Q 2 



G P M {Q 2 ) 



(3.1) 



0.2. 



(3.2) 



The latter and the 



We assume the same C for an estimate of Ff from Ff: Eq. (|3.1|l then provides F^L-, 2 (x,Q 2 ). 
computed f PN < A are from Eq. I|2.4|l seen to be input for the calculation of nuclear SF 

Note that Eq. I|3.2|> for C(l, Q 2 ) requires also G^ 4 , the very FFs we wish to extract. Fortunately, for not too large 
Q 2 the ratio in Eq. I|3.2|) depends only marginally on G 7 ^. Solely in order to determine C(l, Q 2 ), we use in Eq. I|3.2|) 
the parametrization of G 7 ^ given in Ref. l36t we do not demand self-consistency. 

5) Electromagnetic FFs: In Ref. Q we had adopted the ratio 
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G" f 



(G p M fi P ) 
1 

[1 - 0.14(Q 2 - 0.3)] for 0.3 < Q 2 < 5.5 GeV 2 



for Q 2 < 0.3 GeV 2 



choice I 



(3.3) 



from polarization transfer in e(p,p)e |lllll2j|. a p = G P E / (n P G d) ■ 

For all but the smallest Q 2 , the above mentioned E/M ratio for the p FFs's disagrees substantially from results 
from Rosenbluth-separated elastic ep data (Ul^. Two-photon exchange contributions have been computed [T^|). 
but those have as yet not been extracted from data in order to re-analyze the extraction of FFs. 

It is therefore impossible at this moment to make an impartial choice for p FFs, and we shall report results for 
two sets of extracted a n : 'I', based on the ratio (|3.3|l and 'II', using cross section data. We adopt the recommended 
parameterizations of G P E M 



(Q 2 ) 



G p E ,m, 



AQ 



2\ml 



choice II 



(3.4) 
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A previously used parametrization for a p (3^ is close to one of the form (|3.4[1 (see Table 2 in Ref. l36|) . 

There is still lacking reliable data for G 7 ^ beyond relatively low Q 2 [4(| . Analyses continue to prefer a Galster-like 
form 

G n E (Q 2 ) = -HnG d (Q 2 )-^-, (3.5) 
1 + Br] 

with A = 0.942,5 = 4.61 IM El El- Two recent, more precise measurements yield A = 0.888,5 = 3.21 
Differences do no amount too much for the small measured Q 2 : for larger Q 2 , where G% grows in relative importance, 
there is no way to prefer one particular assumed extrapolation. 



IV. GENERAL OBSERVATIONS AND EXTRACTION PROCEDURE. 



We start with the NE parts of inclusive cross sections and remark on two samples for quite different Q 2 , namely 4 He 
(E = 3.6 GeV, (9 = 15°, Q 2 « 0.973 GeV 2 ) and D (E = 4.045 GeV, 6 = 30°, Q 2 w 2.77 GeV 2 ). Figs. 1,2 show on a 
linear scale the cross section data (solid circles), the calculated NI background (dotted curves) and a A ' NE computed 
for some a n , and two additional values, differing from the central one by « 5%. Those curves can be compared with 
the difference between the cross section data and the computed NI component (open circles). The dot-dashed curves 
are empirical NI components, which will be discussed below. The following observations from the above comparison 
are not limited to the chosen examples, but hold quite generally: 

a) Whether a QEP stands out in the data or not, in regions away from x ~ 1, 

a A,NE ig 

barely modified when a n is 

changed by as much as 5% around some average value. Only in the immediate neighborhood of the QEP is <j A ' NE 
sensitive to small variations in a n . 

b) In view of a), a precision extraction of a A ' NE in the region x ~ 1, and thus indirectly of a„, requires a well- 
determined, locally small, inelastic background. The accuracy of such an extraction is limited by the precision in the 
input nucleon SF, the quality of data and of the calculated f PN < A . 

We summarize expectations for inclusive cross sections 0: 

i) On the low-^, elastic side side of the QE region, cross sections for small Q 2 are predominantly NE and the quality 
of the extracted a n depends on the precision with which one can calculate the point-nucleon nuclear SF f PN > A in the 
small wings, away from the peak. 

ii) Approaching the QE region, NE components still dominate, provided Q 2 is not too large. If f PN > A is sharply 
peaked, as is the case for the lightest nuclei, the same will be observed in <r A - ex P w a A ' NE . No matter what the 
kinematic conditions are, the QE region close to x ps 1 is most sensitive to a n and is therefore the primary source for 
the looked-for information. 

hi) On the inelastic, large-;/ side of the QE region, NE components decrease relatively to the increasing inelastic 
background. 

iv) For further increasing v, inclusive excitation of N resonances produces pseudo-resonances in the target. For 
jPN,a w jt n a sharp maximum, pseudo-resonances in nuclei and genuine resonances in nucleons, peak at about the 
same x. For increasing Q 2 , pseudo-resonances gradually coalesce in a smooth background, due to overlapping tails of 
higher pseudo-resonances. 

The extraction procedure follows from the above. First we define reduced cross sections on a p and a composite 
target A 



K v ' (E :9,v) = o M {E- 0, u) 



Fi> A (x,Q 2 ) , 2F?> A (x,Q 2 



M 



-tg>/2) 



(4.1) 



Eqs. (2.7) and (2.8) show that K A ' NE (x, Q 2 ) and f PN < A have a similar rapid variation with x. 

Next we subtract the computed NI background from the data over the entire measured x(z^)-range. We then search 
for a continuous x range for which the difference of reduced cross sections AK A = K A ' exp — K A ' NI and f PN > A shows 
maximal similarity. In that range we identify 

j^A,exp tsA,NI j L tsA,NE 



K J 



and determine for each data-point in the selected x-range 

«n((9 2 )|^n| : 



2MK A > NE /w A v A 



K J 



B 2 h 



\ + tg 2 {6/2)/v A 



2MR A ' NE '- A 



/w 



1/2 



1/2 



[a p (Q 2 )(i p ] 



x,Q 2 



(4.2) 

(4.3) 
(4.4) 
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Above w A (x, Q 2 ) = f PN ' A (x, Q 2 )G 2 (Q 2 ); v A (x, Q 2 ) = x 2 C A (x, Q 2 )/[2(l + V )] and B 2 (Q 2 ) = 2[G£(Q 2 )] 2 /G 2 (Q 2 ). 

We remark that Eq. 12.41) also predicts that K A ' NE and f PN > A have practically the same A-dependence. Eqs. 
(14.3(1 . (|4.4(l show that residual dependence of K A ' NE on both A and cc is far weaker than the same in 

fPN,A In 

addition one notices in Eq. 1)4. 3|) a dependence on 9. For very forward scattering angles the above extraction method 
may occasionally become unstable. An example will be mentioned below. 

Strict fulfillment of the above requirement implies that a n is independent of the data points, which have been 
selected for the extraction. In practice one deals with data for fixed E, 9 and varying x, hence with Q 2 varying 
over the measured ^-range. In general that variation is mild, but not insignificant for Q 2 < 0.7 — 0.8 GeV 2 , and 
the extracted a n will vary there with Q 2 . But even for data sets with fixed Q 2 , experimental inaccuracies and the 
imperfections in the theoretical model, cause extracted a n to depend on the selected data points Xj. Ultimately one 
has for each data set to determine an average (q„). 

It is virtually impossible to incorporate in the analysis all experimental errors and uncertainties in both input and 
theory, and we therefore limit ourselves to the following: 

a) Published tabulated cross sections give statistical and occasionally systematic errors, but frequently only the 
former are shown in figures. Only those are incorporated in our analysis. Most abundant and accurate are D data, 
and occasionally one can assign practically constant relative errors for selected x- intervals. 

b) As mentioned in the Introduction the issue of the proton form factors is n ot y et settled. Two-photon exchange 
contributions to ep inclusive scattering apparently influence on the few % level [Olllj]- Although the E/M ratio for 
the proton as measured in the polarization transfer measurement ls believed to be correct, p FFs, cannot be 
determined, without first to apply the above corrections to cross section data. This has not yet been done. The above 
reflects on both the proton data and the parametrization of Ref. (3^ for a ra , which we used in Eqs. (|4.3[l . (|4.4|) . 

At this point we remark that two-photon exchange contributions cannot be accommodated in a single generalized 
convolution (2.4). The observation, that data for K exp (x) — K NI (x) and f PN ' A (x) have very similar variation with 
x, does not allow more than a few % two-photon exchange contributions. 

We therefore stuck to the procedure followed in Ref. using two sets of proton FFs: one using the Jones results 
and a second one suggested in Ref. [36|. both not yet corrected for two-photon exchange contributions. We expect 
the two methods to provide extremes between which the correct results will ultimately fall. For both sets we applied 
the error analysis a). 

c) The uncertainty in the electric form factor of the neutron appears to be of no consequence. Assuming that the 
different parametrizations may be extrapolated to the largest Q 2 needed, the ratio \G r ^/G p E ] 2 in B 2 , Eq. I|4.3f) may 
vary by as much as 30%. However, for Q 2 > 0.5 GeV 2 the above B 2 term is far smaller than the first term in the 
numerator in Eq. I|4.3(l . and its inclusion affects a n by less than 1%! 

We conclude this Section by the following remark. Our extraction method for (a n ) rests on a test, checking whether 
the difference K A ' NE = K A ' exp — K A ' NI and the SF f PN ' A have similar x-dependence over a continuous set of data 
points. Having determined from those values a n {x{) and their average (a n ), we calculate the corresponding NE 
component NE((a n )) = a A ' NE ({a n }). Likewise Nl(comp) stands for the computed a A ' NI . 

The above defined NE((a„)) are constructed to fit in the mean NE(extr)=data-NI(comp) in the selected range of 
data points Xi around the QEP. Those actually continues to approximately reproduce NE(extr) over a range beyond 
the chosen points of extraction. 

Barring the effect of a mildly varying Q 2 over the the points of a data set, perfect data and an exact theory for NI 
ought to produce a NE((a ra }), fitting NE(extr) over the entire x or v interval. In the following Section we shall find 
that deviations set in beyond some is, where NI about overtakes NE. Those deviations reach a maximum around the 
position of the first pseudo-resonance and then rapidly decreases to 0. The culprit may well be Nl(comp), in which 
case we define an empirical NI by 

Nl(emp) w data - NE((a„)) (4.5) 

By construction NE((a n )) closely fits data-NI(emp). We shall return to a possible source of the apparent insufficiency 
of Nl(comp). 

V. DATA AND RESULTS. 

Below we report on Qf„(G^f), extracted from abundant cross sections for total inclusive scattering of unpolarizcd 
electrons on D and 4 He. Additional information comes from, partly re-analyzed sparse data on their transverse 
components for both targets. 

We start with particulars on data and results collected in Table I. Columns refer to target, beam energy, scattering 
angle, range of x, Q 2 and the value Q 2 = Q 2 {vqep) at the QEP. In the last column we first mention the number of 
selected data points Xj for each data set. Those are followed by the weighted averages of the extracted a n with their 
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error of the mean, which includes statistical errors on the cross sections for both parameterizations I and II discussed 
in Section III. Since systematic errors have been disregarded, the stated error bars should be considered as lower 
limits. 

Only a sample of analyzed data sets are presented in Figs. 1-12. The two options I, II for p FFS produce the same 
elastic ep cross sections with different E/M ratios. For that reason there is no need to specify the option in Figs.: it 
enters only in the ultimately extracted a n - 

HE1): E = 2.02 GeV, 9 = 20°; E = 3.595 GeV, 9 = 16,20°; Q 2 = 0.434,0.873,1.270 GeV 2 0. 

The above NE3 data have previously been analyzed pj and the very good fits on the elastic side of the QEP have 
been attributed to accurate knowledge of the required density matrices and the inclusion of C22 ■ The present refined 
calculations produce the observed rise of the first two cross sections towards the pseudo-resonances. Those are locally 
« 20% short of the first two data sets and « 10% for the third one. Fig. 3 illustrates the latter one. Since the involved 
Q 2 are relatively small, FFs for both options I, II are essentially the same. 

HE2): E = 2.7, 3.3, 3.6, 3.9, 4.3 GeV, 6 = 15°, Q 2 = 0.453, 0.662, 0.781, 0.907, 1.090 GeV 2 ; E = 0.9,1.1 GeV, 
9 = 85°, Q 2 = 0.78, 1.09 GeV 2 . 0,0]. 

The above NE9 cross section data are unpublished parts of the PhD. Thesis of J-P Chen 2]. Those are in principle 
a welcome source of information on a n over a dense Q 2 -range, which partly overlap the Q 2 range of the the NE3 data, 
but extend beyond the QEP and the adjacent minimum, and occasionally into the pseudo-resonance region. 

A comparison with data, illustrated by Fig. 4 for E = 4.3 GeV, 9 — 15° shows a pattern, similar to that for the 
NE3 data. There is a deficiency of « 30% at the peak of the first pseudo- resonance for the set with lowest Q 2 , which 
is reduced to « 10% for the larger Q 2 . 

HE3): Rl.t f° r approximately constant |q| =1.05 and running v @,01- 

Eq. H4.4f) . using Rt is a simpler expression than Eq. (|4.3|) . for reduced, total inclusive cross sections and requires 
only additional knowledge of a p . Extracted a n from, in principle favored Rt data, ought to be close to the ones from 
cross section data HE2) with approximately the same Q 2 , yet Table I shows fairly large deviations. The following 
may well be one of the causes. 

In order to be eligible as partners for a Rosenbluth extraction, some data sets HE2) have been chosen for fixed 
9 = 85° and at two beam energies, such that the x, Q 2 approximately coincide with those of the first set with 9 = 15° 
at different beam energies. Since such a match is never perfect, extrapolations of data are necessary. 

Rt is relatively large for x ~ 1, but that is also region where the point-nucleon SF f PN ' A (x, Q 2 ) varies sharply with 
x. A 3% error or uncertainty in an extrapolation of data points to values for x ~ 1, may cause a 10% change in j PN ' A 
on which F 2 A ' NE depends linearly. We thus tend to actually trust more the involved total cross section information. 
A PWIA analysis of Rt is reported to be in good agreement with data on the elastic side of the QEP ^| . 

The following D data sets are nearly all for appreciably larger Q 2 than for He. 

Dl): E = 4.045 GeV; 9 = 15,23,30,37,45,55°; Q 2 = 0.972, 1.940, 2.774, 3.535, 4.251, 4.900 GeV 2 HEE3 

In our previous analysis we searched for cross section data dominated by Ff . This limited the suitable data to the 
two lowest 9 sets above pj. With reliable parameterizations for both Ff 2 , the above restriction is lifted, providing a 
far larger number of data points for analysis. 

Only in the data for the lowest two angles is a pseudo-resonance clearly visible: For 9 > 30° the pseudo-resonance 
structure gradually disappears. One still notices a minimum beyond the QEP, but for further growing Q 2 , there 
remains no more than a break in the slope of ln( a). The agreement over the entire range of v, with the exception of 
thcintermediate range, discussed at the end of Section IV, is good and sometimes excellent. It comprises qualitatively 
different features, as are the position and intensity at the minimum between the QEP and the pseudo-resonance, as 
well as the position of the peak of the pseudo-resonance. In those regions the NI components first compete with NE, 
overtake those and finally dominate. 

An underestimate of the computed NI in the intermediate ^-range, and its remarkably similar relative size is 
apparent in virtually all analyses. For example, Figs. 5, 2 and 6 for 9 = 15, 30, 55° show a bell-shaped excess over 
NE(extr) on the inelastic side of the QEP, as extracted from Eqs. (|2.1|) . (|2.7|l . (|2.8|l . Dot-dashes show Nl(emp), which 
produce NE((a„)) for all v. 

D2): E = 5.507 GeV, 9 = 15.15,18.98,22.81,26.82°; Q 2 = 1.75,2.50,3.25,4.00 GeV 2 HHJ. 

The kinematics of the above NE11 data partly overlap those of Dl). Virtually all remarks on Dl) hold also for D2) 
(cf. Figs. 7,8 for 9 = 15.15, 26.81°.) 

D3): R T for the x, Q 2 kinematics of D2 [i^. 

For kinematics close to those in D2), the extracted a n from Rt are comparable or slightly higher than from 
D2). Otherwise much the same remarks as for HE3) above hold also here: In order to obtain x,Q 2 matching for a 
Rosenbluth separation, the unavoidable handling of data around the QEP x ~ 1 is bound to be imprecise in view of 
the sensitivity of f PN > D around the QEP, which grows with Q 2 . 

D4): E = 2.015 GeV, 9 = 38.84°; E = 3.8 GeV, 9 = 47.86°; E = 4.212 GeV, 9 = 53.39°; E = 5.12 GeV, 
9 = 56.64°. Q 2 = 1.22,3.15,5.10,6.83 GeV 2 ji| 
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The above NE18 SLAC data for D are for high E and relatively large 9 and are restricted to the immediate QE 
region. In spite of considerable experimental noise, those data are of interest, in view of the large Q 2 involved. 
Computed results are in agreement with data around the QEP, which show much scatter. Beyond that region, 
disagreements are less than rj 25% (cf. Fig. 9). For reasons already mentioned we did not analyze parallel data for 
heavier targets. 

D5): E = 9.744, 12.565, 15.730, 17.301, 18.476, 20.999 GeV, 9 = 10°. Q 2 = 2.5, 4.0, 6.0, 7.1, 8.0, 10.0 GeV 2 [Ell- 
Cross sections computed with rather primitive input show on a tight logarithmic scale reasonable agreement with 
the above old data, but on a linear scale considerable scatter in data is apparent. The above Rock data for small 
9 provide a unique example of marginal stability: the E — 9.744 provide two adjacent subsets with rather different 
average for a n . Both cause (a n ) to have relatively large error bars (Table I). 

The main interest is their values out to the largest Q 2 , measured until this date. Figs. 10-12 show results for 
E = 15.73, 18.476, 20.999 GeV. Table I shows that for E = 18.476 (a n ) is relatively large, and not in line with other 
(a n ) in the data set. One observes close correspondence between (a n ), extracted from different data sets. 

The overall outcome is compiled in Figs. 13-15. In inserts we give symbols for experiments (empty ones for older 
and filled symbols for re-analyzed ones), kind if data and Ref. numbers. Fig. 13 contains results from He data, which 
are all for low Q 2 < 1.3 GeV 2 with a n , close to 1. For those Q 2 , differences for the choices I, II are relatively small 
(see Table I); displayed results are for I. In detail: 

i) Asymmetry measurements in QE 3 He(e, e') |5lll53. IH^ . 

ii) Preliminary results of the analysis of i) , without 3-body channels in the FSI Ej3 . 

iii) Ratios of exclusive D break-up cross sections, with p, respectively n detected |55l l56l l57| (In view of the criticism 
of Jourdan et al. |5^. we omit some old data |59jl. 

iv) The above HE1)-HE3) HUB 3 He(e,e'). 

Fig. 14 contains results, extracted from the sets Dl)-D5) with Q 2 > 1.0 GeV 2 H E EHH |H H3 ■ The only 
displayed result of a previous analysis is D2) |4^. Fig. 15 is the same as Fig. 14 for choice II of FFs. 
In both the drawn lines are inverse polynomial fits of a n (Q 2 ) from the D data, as used in Ref. j3|| 

a n (Q 2 f = l/[l + 0.007437Q 2 + 0.002815161Q 4 - 0.000115008O 6 ] 

a n {Q 2 ) H = 1/[1 + 0.06568Q 2 - 0.00678662Q 4 + 0.000323355Q 6 ] (5.1) 

Although only approximately valid for Q 2 > 1.5 GeV 2 ), the fits extrapolate to the 'correct' a n (0) = 1. In spite of 
scatter in the data two observations on (a n ) stand out: 

a) a n (Q 2 ) clearly decreases with increasing Q 2 

b) There is close correspondence between extracted values from different data sets, for instance for Q 2 = 2.5,4.0, ~ 
7.0 GeV 2 . 



VI. SUMMARY AND CONCLUSIONS. 



The overwhelming majority of its properties cannot been measured on free neutrons and one is therefore led to 
study those on neutrons, which are bound in nuclear targets. The extraction of any such property thus demands an 
accurate treatment of the embedding of the neutron in that target. Even when feasible, one has in addition to assume, 
that those extracted quantities are the same as for a neutron in vacuum, i.e. that the former are not intrinsically 
modified by the medium. Nowadays one can accurately compute nuclear properties for A < 4, and there seems to be 
no evidence for medium effects. Those may well be artifacts of approximations. 

The present contribution deals with the extraction of the reduced magnetic form factor of the neutron from cross 
sections for the inclusive scattering of electrons from the lightest targets. 

Part of the available data have been analyzed before [lj . There we limited ourselves to those kinematic parts of 
data sets, where the elastic absorption of an exchanged photon on a nucleon had been estimated to be much in excess 
of the inelastic ones, which were subsequently disregarded: That procedure not only affected accuracy. It also limited 
the analysis to the immediate neighborhood of QE peaks, where the NI background is very small. That region covers 
only a small section of the data, which generally stretch over wide energy-loss ranges. 

Ref. followed the principal ideas of Lung and of older work, using rather primitive tools and input jE^j. 
A number of incentives called for a re-analysis of the above results. Together with our wish to include the above 
mentioned neglected kinematic areas, we included previously disregarded data. We used moreover recently published, 
precise parameterizations of the input SF Ff 2 for the proton, valid through resonances and reaching into the DIS 
region. For medium and large Q 2 those parametrizations become unreliable and we had to fall back on parameter- 
izations of resonance-averaged F%. In addition, it appeared possible to up-date the input for the determination of 
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the additionally required neutron SF. Finally we could also sharpen some theoretical tools. As a result the present 
analysis constitutes a significant refinement of the previous one. 

The above comprehensive theory for the extraction of the magnetic form factor of the neutron G^j addressed 
relatively abundant total cross for inclusive electron scattering, and scarce data on their transverse components. In 
principle several targets are accessible to an analysis, but only for the lightest nuclei can one presently calculate with 
great precision nuclear information, which describes the above embedding. Most of the available data on D are of 
good quality and contain several data sets covering a range of partly overlapping Q 2 . The latter fact enables desirable 
consistency checks. For 4 He there are only available rather old, low-Q 2 data of lesser quality and for which also the 
theory is less accurate than for higher Q 2 . Nevertheless we analyzed all. 

The cornerstone of our analysis is the possibility to reliably compute components of the nuclear inclusive cross 
sections, due to inelastic virtual photon absorption on nucleons, provided are known. Those nuclear NI processes 
dominate virtually all kinematic regions, except the ones around the QE peak, where elastic absorption of virtual 
photons on nucleons competes with inelastic processes. 

Cross sections in those QE regions contain the desired information on form factors, and their isolation is therefore 
of primary importance. Simple theoretical considerations predicts the same ^-dependence of the NE components of 
cross sections and the calculable Structure Function f PN ' A (x,Q 2 ) of a fictitious target, composed of point- nucleons. 
The latter drops sharply from its QE peak at x w 1: in the case of a D by a factor w 10 — 50. 

One thus compares the functional dependence on x of the above difference with the same for f PN ' A (x,Q 2 ). In 
regions where close similarity is found, one identifies that difference with the desired NE components of cross sections. 

For all data sets, the above differences between measured reduced total cross sections and calculable inelastic 
backgrounds, appear to follow the predicted x-behaviour, roughly for x > 0.7—0.8. The above provides incontrovertible 
proof that in the above restricted areas, where the x-dependence is most outspoken, the above defined differences are 
indeed the NE components of the cross sections. 

For each data set we then selected a continuous x-range, for which the correspondence with the x-dependence of 
jPN,a j s b es t_ Dependent on the quality of the data set, the number of thus selected x points may be as large as 17. 
Ideally, those should reproduce G^, independent of x or A. 

We first summarize our results: 

A) There is similarity and sometimes reasonable correspondence between the extracted a n from different targets. 
Differences may in part be due to the same in the quality of the older low Q 2 , He data and the more recent D 
experiments. It is in particular satisfactory that the high-quality E — 4.045 and E — 5.507 D data yield corresponding 
&n(Q 2 ) for similar Q 2 . The same is the case, if the old data of Rock etal. are included. 

B) Our results confirm and reinforce an important conclusion already reached in our previous analysis: a n from 
D data seems to fall faster with Q 2 than the older Lung results indicate. The same is the case for the slope of the 
older He data, but as stated above, we have reservations regarding the used parametrization for F^ for the lowest Q 2 
and the quality of the data. For growing Q 2 , the NE3 and NE9 4 He data and the average behavior of the D data for 
corresponding Q 2 produce similar a n . 

We have emphasized the identification of <j A > NE over finite x-ranges. On the inelastic side of the QEP for virtually all 
D data sets, the use of NI, computed with the recent parametrization of -Flf, produces deviations from NE components, 
computed with (a n ). Those roughly start where NEs=sNI, and grow to a 10-15 % under-estimate towards the peak of 
the pseudo-resonance position, beyond which the discrepancy rapidly disappears. An increase of the NI background 
of just the above size extends the above local fit of the computed NE cross section to comprise the entire z/-range. 

The discrepancy may be the result of a relatively modest under-estimate of the transition strength for inclusive 
excitation of the first resonance for larger Q 2 . It would effect mostly the low v tail of NI, because more inelastic parts 
are screened by overlapping resonances. We emphasized that one cannot apply an ad hoc change of one parameter in 
a multi-parameter fit of F%. 

In spite of the manifestly succesful local isolation of elastic components of the inclusive cross sections, the actual 
extraction of G\ from those is not straightforward. We mention a few sources which may influence the precision of 
the extraction: 

1) Disregarding experimental accuracies, the test of the central requirement is optimal when the inelastic background 
is a small fraction of the total cross section. The latter is the case in the immediate neighborhood of the QEP, where 
the elastic components are most sensitive to variations of G\ ;i . The small values of the inelastic background there 
are due to the tails of overlapping resonances, which by nature have least accuracy. Consequently, changes in a very 
small background around the QEP may have 1% effects on the extracted a n . 

2) Total inclusive cross sections contain in principle Meson Exchange Contributions (MEC), which should be 
eliminated from the data, before those are manipulated as described. Little is known of MEC contributions for high- 
Q 2 inclusive processes. As a measure of their size we suggested recent results for the exclusive processes 3 He(e,e'p)D 
and 4 He(e,e'p) 3 He. Polarization variables appear sizeably affected by MEC, but only tiny corrections are reported on 
cross sections for inclusive scattering with unpolarized beams |6l) . 
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3) Present uncertainties regarding the proton form factors. 

The above concludes our program to extract from data on inclusive cross sections. We are well aware that our 
method is an indirect one, which forced us to make careful checks on intermediate steps. Those unfortunately could 
not circumvent present uncertainties in input. At least the unknown behavior of for large Q 2 seems not to be of 
relevance. Unless the unknown G 7 ^ for larger Q 2 will turn out to deviate strongly from assumed extrapolations from 
low Q 2 , its influence will remain negligible. 

The above clearly demonstrates that the described extraction method is a realistic and internally consistent one. We 
do not think that theoretical tools can be much improved, but it is highly desirable to eliminate uncertainties in some 
input elements of the calculations. In parallel new inclusive scattering data on 3 He and 4 He could add information 
and furnish further proof of consistency. 

In the meantime we are looking forward to final results for G^j from the D(e, e'p)X / D(e, e'n)X experiment of 
Brooks etal. Preliminary results seem to yield a smaller slope for a n than in Figs. 14,15 62] . 
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FIG. 1: Inclusive cross section on 4 He for E = 3.6 GeV, 6 = 15°, Q 2 = 0.781 GeV 2 . Filled and empty circles represent SLAC 
NE9 data Q with their statistical errors, respectively cross sections, from which the theoretical NI part (dotted curve) has 
been subtracted. Empty circles should carry error bars coming from the data and the NI part, but since we cannot at present 
estimate those,we have preferred not show them, up to the last extraction point used. The difference is compared with a Hc ' NE 
for otn = 1.00, 1.00 ± 0.05 (curves). Dot- dashes represent an Nl(emp), defined as a tot — a NE . 
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FIG. 2: As Fig. 1 for D at E = 4.045GeV,6> = 30°, Q 2 = 2.774 GeV 2 . Data: Jlab E89-009 46]. NE curves are for 
a n = 0.95, 0.95 ±0.05. 
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FIG. 3: As Fig. 1 for 4 He at E = 3.595 GeV, 9 = 16°. Q 2 = 0.873 GeV 2 . Data: SLAC- Virginia NE3 |. NE curve for 
a n = 1.006. 
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FIG. 4: As Fig. 1 for 4 He at E = 3.9 GeV, 6 = 15°, Q 2 = 0.907 GeV 2 . Data: SLAC NE9 0. NE curve for a n = 1.007. 
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FIG. 6: As Fig. 1 for D at E = 4.045 GeV, 6 = 55°, Q 2 = 4.900 GeV 2 . Data: JLab [«J. NE curve for a n = 0.991. 
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FIG. 7: As Fig. 1 for D at E = 5.507 GeV, 9 = 15.15°, Q 2 = 2.750 GeV 2 . Data: SLAC NE11 @. NE curve for a n = 1.008. 
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FIG. 8: As Fig. 1 for D at E = 5.507 GeV, 9 = 26.86°, Q 2 = 4.000 GeV 2 . Data: SLAC NE11 @. NE curve for a n = 0.951. 
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FIG. 9: As Fig. 1 for D at E = 5.12 GeV, 9 = 56.64°, Q 2 = 6.83 GeV 2 . Data: SLAC NE18 47]. NE curve for a n = 0.842. 
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FIG. 10: As Fig. 1 for D at E = 15.73 GeV, 9 = 10°, Q 2 = 6.00 GeV 2 . SLAC E133 data NE curve for a n = 0.854 
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FIG. 11: As Fig. 1 for D at E = 17.301 GeV, 6 = 10°. Q 2 = 7.12 GeV 2 . SLAC E133 data [H^- NE curve for a n = 0.879 
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FIG. 12: As Fig. 1 for D at E = 20.999 GeV, = 10°. Q 2 = 10.06 GeV 2 . SLAC E133 data |S(|. NE curve for q„ = 0.769 
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TABLE I: Extraction of a n (Q 2 ) from QE inclusive scattering data on D, 4 He. Columns 1-3 give the target, beam energy E 
and scattering angle 9. Columns 4-6 indicate the ranges of u,x,Q 2 , selected for the analysis of the QE and pseudo-resonance 
region, and Q 2 at the QEP. The last column gives the number of points, selected for the extraction of (a n (Q 2 )), its average 
(with Q 2 — > Q 2 for identification and error of the mean for the new, respectively the old FF parametrization. 



target 


E [GeV] 


6 


v[GeV] 


X 


Q 2 [GeV 2 ]; Q 2 


n [(a n (Q 2 ))o]\ [K«? 2 ))iv]] 


4 Hc [4] 


2.02 
3.595 
3.595 


20° 
16° 
20° 


0.120-0.435 
0.180-0.750 
0.255-0.870 


2.054-0.473 
2.814-0.563 
3.024-0.723 


0.463-0.386; 0.434 
0.951-0.792; 0.873 
1.450-1.180; 1.270 


5 1.011±0.032; 1.008±0.032 
11 1.006±0.045; 0.954±0.045 
8 0.940±0.031; 0.865±0.030 


4 Hc [2] 


2.7 
3.3 
3.6 
3.9 
4.3 


15° 
15° 
15° 
15° 
15° 


0.139-0.784 
0.099-0.858 
0.100-0.793 
0.104-0.698 
0.191-1.000 


1.806-0.239 
3.872-0.341 
4.573-0.462 
5.167-0.649 
3.357-0.515 


0.471-0.352; 0.453 
0.720-0.549; 0.662 
0.859-0.689; 0.781 
1.009-0.851; 0.907 
1.204-0.967; 1.090 


4 1.120±0.035; 1.114±0.036 
4 1.064±0.029; 1.031±0.029 
7 1.028±0.037; 0.982±0.038 
10 1.007±0.034; 0.950±0.034 
9 0.978±0.033; 0.912±0.033 


4 Hc [2] 

mm! 


0.9 
1.1 


85° 
85° 


0.316-0.696 
0.492-0.796 


1.617-0.256 
1.322-0.408 


0.959-0.335; 0.78 
1.221-0.610; 1.09 


5 1.014±0.030; 1.009±0.029 
5 0.994±0.062; 0.979±0.062 


DHc 






0.390-0.600 


1.298-0.659 


0.950-0.741; 0.90 


5 1.079±0.006; 1.066±0.006 


£[2, 3] 


4.045 
4.045 
4.045 
4.045 
4.045 
4.045 


15° 
23° 
30° 
37° 
45° 
55° 


0.405-1.185 
0.615-1.575 
0.930-2.505 
1.245-1.935 
1.695-2.505 
2.160-3.045 


1.320-0.355 
1.911-0.537 
1.933-0.355 
1.951-0.946 
1.750-0.776 
1.603-0.603 


1.003-0.789; 0.972 
2.206-1.589; 1.940 
3.376-1.669; 2.774 
4.561-3.437; 3.535 
5.568-3.649; 4.251 
6.502-3.449; 4.900 


10 1.004±0.023; 0.945±0.024 
9 0.994±0.019; 0.921±0.021 
9 0.991±0.017; 0.915±0.019 
9 0.959±0.019; 0.874±0.020 
12 0.936±0.030; 0.845±0.035 
7 0.879±0.018; 0.789±0.018 


D[48] 


5.507 
5.507 
5.507 
5.507 


15.15° 
18.98° 
22.81° 
26.82° 


0.698-1.309 
1.011-1.665 
1.571-2.053 
1.983-2.424 


1.404-0.654 
1.418-0.736 
1.149-0.771 
1.121-0.803 


1.840-1.606; 1.750 
2.692-2.300; 2.500 
3.388-2.973; 3.250 
4.176-3.653; 4.000 


11 1.008±0.027; 0.933±0.026 
17 0.998±0.039; 0.919±0.040 
17 0.987±0.049; 0.900±0.051 
6 0.951±0.058; 0.851±0.059 


Rx 






0.773-0.971 
1.281-1.394 
1.692-1.838 
2.095-2.210 


1.066-0.950 
1.051-0.940 
1.034-0.942 
1.027-0.941 


1.769-1.733; 1.750 
2.529-2.461; 2.500 
3.286-3.188; 3.250 
4.059-3.905; 4.000 


10 0.957±0.023; 0.940±0.024 
5 0.973±0.022; 0.957±0.022 
7 0.978±0.018; 0.959±0.018 
4 0.901±0.025; 0.866±0.025 


£>[4JJ 


2.015 
3.188 
4.212 
5.120 


38.84° 
47.68° 
53.39° 
56.64° 


0.582-0.722 
1.635-1.786 
2.659-2.810 
3.567-3.719 


1.168-0.850 
1.054-0.871 
1.058-0.904 
1.069-0.928 


1.277-1.152; 1.22 
3.235-2.921; 3.15 
5.280-4.767; 5.10 
7.156-6.461; 6.83 


7 0.924±0.033; 0.868±0.036 
4 0.956±0.089; 0.897±0.096 
6 0.855±0.064; 0.774±0.070 
4 0.842±0.051; 0.786±0.049 




9.744 
12.565 
15.730 
17.301 
18.476 
20.999 


10° 
10° 
10° 
10° 
10° 
10° 


1.14-1.60 
1.821-2.537 
2.758-3.768 
3.274-4.331 
3.087-4.748 
4.658-5.885 


1.19-0.803 
1.200-0.804 
1.197-0.808 
1.200-0.838 

1.191- 0.864 

1.192- 0.873 


2.547-2.4112; 2.50 
4.101-3.828; 4.00 
6.199-5.716; 6.00 
7.373-6.817; 7.12 
8.289-7.705; 8.03 
10.424-9.642; 10.06 


6 0.977±0.056; 0.878±0.061 

7 0.913±0.050; 0.784±0.053 
10 0.854±0.029; 0.714±0.031 

8 0.879±0.124; 0.812±0.117 

9 0.958±0.038; 0.912±0.038 
5 0.769±0.118; 0.708±0.131 
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FIG. 13: a n (Q 2 ), extracted from recent reliable He data, using choice I, (3.3) for G P E . See last paragraph of Section V for 
description. 
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FIG. 14: Same as Fig. 14 from D data data for choice I. 
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FIG. 15: Same as Fig. 14, using choice II (3.4) for G" E . 



